Central limits and homogenization in random media 



Guillaume Bal * 
February 2, 2008 

Abstract 

We consider the perturbation of elliptic operators of the form P(x, D) by 
random, rapidly varying, sufficiently mixing, potentials of the form q{j,ij). We 
analyze the source and spectral problems associated to such operators and show 
that the properly renormalized difference between the perturbed and unperturbed 
solutions may be written asymptotically as e — > as explicit Gaussian processes. 
Such results may be seen as central limit corrections to the homogenization (law 
of large numbers) process. Similar results are derived for more general elliptic 
equations in one dimension of space. 

The results are based on the availability of a rapidly converging integral formu- 
lation for the perturbed solutions and on the use of classical central limit results 
for random processes with appropriate mixing conditions. 

keywords: Homogenization, central limit, mixing coefficients, partial differential equa- 
tions with random coefficients, random oscillatory integrals. 

AMS: 35R60, 35J05, 35P20, 60H05. 

1 Introduction 

There are many practical applications of partial differential equations with coefficients 
that oscillate at a faster scale than the scale of the domain on which the equation is 
solved. In such settings, it is often necessary to model the rapidly oscillatory coefficients 
as random processes, whose properties are known only at a statistical level. The nu- 
merical simulation of the resulting partial differential equation with random coefficients 
becomes a daunting task. 

Two simplifications are then typically considered. The first simplification consists 
in assuming that the coefficients oscillate very rapidly and replacing the equation with 
random coefficients by a homogenized equation with deterministic (effective medium) 
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coefficients. The homogenization of many linear and nonlinear partial differential equa- 
tions with periodic [71 [39] and random highly oscillatory coefficients has been obtained 
to date [m [la [291 [Ml [Ml [371 ig. 

The solution to the equation with random equations may also be interpreted as a 
functional of an infinite number of random variables and expanded in polynomial chaoses 
[T6l W7\ . A second simplification consists then in discretizing the randomness in the co- 
efficients over sufficiently low dimensional subspaces -primarily by Galerkin projection- 
so that the partial differential equation with random coefficients may be fully solved 
numerically. We refer the reader to e.g. [31 [221 [2S1 [3S1 US] for references on this 
active area of research. Such problems, which are posed in domains of dimension d + Q, 
where d is spatial dimension and Q the dimension of the random space, are computa- 
tionally very intensive, although they have the main advantage of providing realistic 
statistical fluctuations of the random solution, which are absent in the homogenization 
approximation. 

The two aforementioned approaches can hardly be reconciled. Homogenization arises 
in a limit where the law of large numbers applies and the solution becomes asymptoti- 
cally a deterministic quantity. The number of random variables describing the random 
coefficients thus tends to infinity, a limit that is difficult to obtain by polynomial chaos- 
type expansions. 

In several practical settings such as e.g. the analysis of geological basins or the 
manufacturing of composite materials, one may be interested in an intermediate sit- 
uation. We may observe experimental fluctuations in the random solution which are 
not accounted for by the homogenized solution, and yet, may be in the presence of a 
sufficiently rich random environment so that full solutions of the equation with random 
coefficients may not be feasible. This is the type of settings that motivate the studies 
of this paper. Our main objective will be to characterize the statistical structure of 
the corrector to the homogenized, deterministic, limit. Whereas the deterministic limit 
may be seen as a law of large number effect, we are interested in characterizing the next 
order term, which arises as an application of the central limit theory. 

In most practical cases of interest, starting with the elliptic operator V ■ ae(x, ct;)V, 
with X G -D C M'^ and u ^ Q the space of random realizations, the calculation of 
the homogenized tensor is difficult and does not admit analytic expressions except in 
very simple cases [29]. The amplitude of the corrector to homogenization, let alone 
its statistical description, remains largely open. The best estimates currently available 
in spatial dimension d > 2 may be found in [IS]; see also [TSl [20], [H] for discrete 
equations, and [1] for applications of such error estimates. Only in one dimension of 
space do we have an explicit characterization of the effective diffusion coefficient and of 
the corrector [13]. Unlike the case of periodic media, where the corrector is proportional 
to the size of the cell of periodicity e, the random corrector to the homogenized solution 
is an explicitly characterized Gaussian process of order y/e when the random coefficient 
has integrable correlation [13]. In the case of correlations that are non integrable and 
of the form R{t) ~ for some < a < 1, the corrector may be shown to be still an 
explicitly characterized Gaussian process, but now of order [5]. 

The reason why explicit characterizations of the correctors may be obtained in [5l [13] 
is that the solution to the heterogeneous elliptic equation may be written explicitly. 
Correctors to homogenization have been obtained in more general settings. The analysis 
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of homogenized solutions and central limit correctors to evolution equations with time 
dependent randomly varying coefficients is well known; see e.g. [9], [221 IMl Ell [351 HI] 
for reference on the Markov diffusion approximation and the method of the perturbed 
test function. In the context of the one-dimensional Helmholtz equation, this would 
correspond to solving the equation on an interval (0, a) with initial conditions of the 
form Ue{0) and ^^(O) known. The asymptotic limit of boundary value problems, which 
corresponds in our example to prescribing Me(0) and Me (a), requires somewhat different 
mathematical techniques. We refer the reader to [251 US] for results in the setting 
of one-dimensional problems. Note that in the case of a much stronger potential, in 
dimension d = 1 of the form e'^q^ instead of in the above Helmholtz operator, the 
deterministic homogenization limit no longer holds. Somewhat surprisingly, the solution 
of a corresponding evolution equation still converges to a well identified limit; see |43j . 

In spatial dimensions two and higher, a methodology to compute the Gaussian fluctu- 
ations for boundary value problems of the form —Aue+F{u£, x, -) = /(x) was developed 
in [23] . An explicit expression for the fluctuations was obtained and proved for the linear 
equation (— A + A + g(-))Me(x) = /(x) in dimension d = 3. In this paper, we revisit the 
problem and generalize it to linear problems of the form P(x, D)^^ + ^^(x)^^ = /(x) 
with an unperturbed equation P(x, D)m = /, which admits a Green's function G(x, y) 
that is more than square integrable (see below). The prototypical example of interest 
is the operator P(x, D) = —V ■ a(x) V + q'o(x) with sufficiently smooth (deterministic) 
coefficients a(x) and qo{'x) posed on a bounded domain with, say, Dirichlet boundary 
conditions, for which the Green's function is more than square integrable in dimensions 
l<d<3. 

Under appropriate mixing conditions on the random process ^^(x, u;), we will show 
that arbitrary spatial moments of the correctors 



where and u are the solutions to perturbed and unperturbed equations, respectively, 
and where M is a smooth function, converge in distribution to Gaussian random vari- 
ables, which admit a convenient and explicit representation as a stochastic integral with 
respect to a standard (multi-parameter) Wiener process. If we denote by ui the weak 
limit of = e~2[ui, — uq), we observe, for 1 < d < 3, that E{f ^^(x, cj)} converges 
to E{m^(x, a;)}, where is the leading term in up to an error term we prove is 
of order 0{e'^) in L^{fl x D). This shows that the limiting process ui captures all the 
fluctuations of the corrector to homogenization. This result is in sharp contract to the 
cases d > 4 and to homogenization in periodic media in arbitrary dimension, where the 
weak limit of the corrector captures a fraction of the energy of that corrector. We thus 
see that corrections to homogenization are somewhat different for 1 < d < 3 and d > 4 
so that the square integrability of the Green's function is a natural condition in the 
framework of homogenization in random media. 

We obtain similar expressions for the spectral elements of the perturbed elliptic 
equation. We find that the correctors to the eigenvalues and the spatial moments of 
the correctors to the corresponding eigenvectors converge in distribution to Gaussian 
variables as the correlation length e vanishes. In the setting d = 1, we obtain similar 
result for more general elliptic operators of the form — T"(ae j-) + Qo + qe by appropriate 
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use of harmonic coordinates [3l] . The extension to similar operators in dimension d > 2 
remains open. 

As was already mentioned, the theory developed here allows us to characterize the 
statistical properties of the solutions to equations with random coefficients in the limit 
where the correlation length (the scale of the heterogeneities) is small compared to the 
overall size of the domain. In many practical problems, it is a good approximation to 
the statistical structure of the solution of the equation and possibly be all that one is 
interested in. 

Asymptotically explicit expressions for the correctors may also find applications in 
the testing of numerical algorithms. Several numerical schemes have been developed 
to estimate the heterogeneous solution accurately in the regime of validity of homog- 
enization by using discretizations with a length scale h that is large compared to the 
correlation length: h ^ S] see e.g. [H [21 [20l [21], [IQ]. A possible application of the 
explicit expression for the correctors is to see whether these algorithms can capture the 
central limit correction to the solutions to the random partial differential equations. 

An other application concerns the reconstruction of the constitutive parameters of 
a differential equation from various measurements, for instance the reconstruction of 
the potential in a Helmholtz equation from spectral measurements [30l [H]. In such 
cases, reconstruction algorithms provide lower-variance reconstructions when the cross- 
correlations are known and used optimally in the inversion; see e.g. Provided that 
go is the deterministic quantity that we wish to reconstruct, higher-frequency compo- 
nents that we may not hope to reconstruct still influence available data. The correctors 
obtained in this paper provide asymptotic estimates for the cross-correlation of the 
measured data, which allow us to obtain lower- variance reconstructions for go! see [6]. 

An outline for the rest of the paper is as follows. Section [2] considers the convergence 
of the corrector to the homogenized solution for the Helmholtz equation with source 
term in dimensions 1 < ci < 3. The proof is based on showing the rapid convergence of 
a Lippman-Schwinger-type integral formulation (see ([8]) below), and on applying central 
limit theorems to random oscillatory integrals. The behavior of the oscillatory integrals 
is considered in arbitrary dimensions in section [231 where a comparison between homog- 
enization in random and periodic environments is also considered. The generalization 
to a more general one-dimensional elliptic source problem in detailed in section [31 The 
results on the correctors obtained for source problems are then extended to correctors 
for spectral problems in section [H The proof is based on adapting classical results [31] 
on the convergence of the spectra of operators that converge on average in the uniform 
norm. The results obtained for the spectral problems are then briefly applied to the 
analysis of evolution equations. Some concluding remarks are presented in section [3 

2 Correctors for Helmholtz equations 

Consider an equation of the form: 




(1) 
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where P(x, D) is a (deterministic) self-adjoint, elliptic, pseudo-differential operator and 
D an open bounded domain in R*^. We assume that P(x, D) is invertible with symmetric 
and "more than square integrable" Green's function. More precisely, we assume that 
the equation 

P(x,D)n = /, xgD 

M = X G ^ ' 

admits a unique solution 

u(x) = ^/(x):= / G(x,y)/(y)dy, (3) 

JD 

and that the real-valued and non-negative (to simplify notation) symmetric kernel 
^(x, y) = G{y, x) has more than square integrable singularities so that 

X 1-^ ^ y [(^^"^"^(x, y)(iyj is bounded on D for some r/ > 0. (4) 

The assumption is typically satisfied for operators of the form P(x, D) = —V ■ a(x) V + 
(t(x) for a(x) uniformly bounded and coercive, o"(x) > 0, and in dimension d < 3, with 
r] = +00 when d = 1 (i.e., the Green's function is bounded), 77 < cxd for (i = 2, and 
rj < 1 for d = 3. 

Let ge(x, cj) = q{-,uj) be a mean zero, (strictly) stationary, process defined on an 
abstract probability space (f2,jF, P) [TJ]. The process ^^(x, cj) will be modified in the 
sequel as the process ge(x, u) appearing in ([1]) to ensure that solutions to the Helmholtz 
equation exist. We assume that g(x, cj) has an integrable correlation function: 

P(x) =E{g(y,cu)g(y + x,a;)}, (5) 

where E is mathematical expectation associated to P. The above expression is indepen- 
dent of y by stationarity of the process g(x, a;). We also assume that g(x, tu) is strongly 
mixing in the following sense. For two Borel sets A,Bc M"', we denote by J-'a and jFg 
the sub-o" algebras of generated by the field g(x, cu). Then we assume the existence 
of a (p— ) mixing coefficient (p{r) such that 



E{(r/-EM)(e-E{0)} 



{E{r^^}E{e}y 



<cp{2d{A,B)) (6) 



for all (real-valued) random variables rj on (fi, JF4, P) and ^ on (f2, jFg, P). Here, d{A, B) 
is the Euclidean distance between the Borel sets A and B. The multiplicative factor 
2 in ([2]) is here only for convenience. Moreover, we assume that (f{r) is bounded and 
decreasing. We will impose additional restrictions on the process to ensure that the 
equation ([1]) admits a solution. 
We formally recast ([1]) as 

Us = G{f - qeUe), (7) 
where Q = P(x, D)~^, and after one more iteration as 

Ue = Gf - QqsQf + GqeGqeUe- (8) 

This is the integral equation we aim to analyze. 
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2.1 Existence and error estimates 

In order for the above equation to admit a unique solution, we need to ensure that 
(/ — QqeGqe) is invertible P— a.s. We modify the process ge(x, u) defined above on a set 
of measure of order e so that Qq^Qq^ has spectral radius bounded by p < 1 P— a.s. To 
do so and to estimate the source terms in ([8]), we need a few lemmas. 

Lemma 2.1 Let q{:x.,uj) he strongly mixing so that ([6]) holds and such that'K{q^} < oo. 
Then, we have: 

|E{g(xi)g(x2)g(x3)g(x4)}| < sup ^H\yi - y3\)^H\y2 - yMq^f' ■ (9) 

{yfc}i<fe<4={xfe}i<fc<4 

Here, we use the notation a <b when there is a positive constant C such that a < Cb. 

Proof. Let yi and y2 be two points in {xfc}i<fc<4 such that c?(yi,y2) > (i(xj,Xj) for 
all 1 < < 4 and such that (i(yi, {zg, Z4}) < rf(y2, {zg, Z4}), where {yi, y2, Z3, Z4} = 

{Xfc}l<fc<4- 

Let us call y^ a point in {z3,Z4} closest to yi. We call y4 the remaining point in 
{x;,.}i<fc<4. We have, using (jS]) and Ejg} = 0, that: 

|E{g(x0g(x2)g(x3)g(x4)}| <^(2|yi-y3|)(E{g2})i(E{(g(y2)g(y3)g(y4))'})i 

The last two terms are bounded by E{g^}^ and E{g®}^, respectively, using Holder's 
inequality. Because (p{r) is assumed to be decreasing, we deduce that 

|E{g(xi)g(x2)g(x3)g(x4)}| < v^dyi - ysD^iq^jl (10) 
If y4 is (one of) the closest point (s) to y2, then the same arguments show that 

|E{g(xi)g(x2)g(x3)g(x4)}| < ^(|y2 - y4|)E{g6}i (11) 
Otherwise, y3 is the closest point to y2, and we find that 

|E{g(xi)g(x2)g(x3)g(x4)}| < ^(2|y2 - y3\)E{q^}l 

However, by construction, |y2 - y4| < |yi - y2| < |yi - yal + Iys - y2| < 2|y2 - ysl, so 
(ITTil is still valid (this is the only place where the factor 2 in ([6]) is used). 

Combining ffTOj) and ffTTj) . the result follows from a At < {ab)^ for a,b > 0, where 
a Ab = min(a, b). □ 

Lemma 2.2 Let q^ be a stationary process ^^(x, u) = g(^, 00) with integrable correlation 
function in Let f be a deterministic square integrable function on D. Then we have: 

mQQeQf\\UD)}<^WLHDy (12) 
Let q^ satisfy one of the following additional hypotheses: 
[HI] g(x, cj) is uniformly bounded F-a.s. 

[H2] E{g^} < 00 and g(x, tu) is strongly mixing with mixing coefficient in ([6]) such that 
(p2[r) is bounded and r'^'^ if 2 (^r) is integrable on M"*". 
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Then we find that 

mG'leQqs\\liLHD))}<^'- (13) 

Remark 2.3 Note the assumption [H2] combined with (p{r) decreasing together impose 
that (f{r) = o(r~^'^). For otherwise, we would have an increasing sequence r„ — ^ oo as 
n — > oo such that V52(r„) > Cr~'^ for some C > 0, and then, since ip2 is also decreasing. 



X r^+i c^r^+i rfr„+i 



Now if there is an infinite number of terms n such that r„+i > 2r„, then there is an 
infinite number of terms such that > ^ and the above sum is infinite. If there 

c<r„+i — 2d 

is a finite number of such terms, then for all n > uq for no sufficiently large, we have 
Tn+i < 2r„ so that 



" n>no n>no "0 



Our assumptions then impose that i^ir) decay faster than r 
Proof [Lemma [2^2] . We denote || ■ || = || ■ \\l'^{d) and calculate 

gqeGfi^) = (^G(x,y)g,(y)G(y,z)t/y)/(z)rfz, 

so that by the Cauchy-Schwarz inequality, we have 

\QqeQf{^)? < ll/ir^ (^G(x,y)g,(y)G(y,z)rfy)'dz. 

By definition of the correlation function, we thus find that 

mOQsGff} < ll/f ^^G(x,y)G(x,C)i?(^)G(y,z)G(C,z)dxrfyrfCrfz. (14) 

Extending G(x, y) by outside D x D, we find in the Fourier domain that 

mGQeQm<\\fr [ I iG(x'rjG(z,-)p(p)5'^%p)rfprfxrfz. 



D2 



Here /(^) = j^^i e *^ ''/(x)(ix is the Fourier transform of /(x). Since -R(x) is integrable, 
then R{e'p) (which is always non-negative by e.g. Bochner's theorem) is bounded by a 
constant we call Rq so that 

mG<leQm < l|/HVi?o / G^(x,y)G2(z,y)rfxrfyrfz < l|/||Vi?o, 

by the square-integrability assumption on G{'x, y). This yields (1121) . Let us now consider 
(ITSi) . We denote by H^ge^gell the norm \\GqeG(le\\c{L'2(D)) and calculate that 

QqeGqe(p{^) = J (yj G'(x,y)ge(y)G(y,z)rfyjge(z)0(z)rfz. 



Therefore, 

(gqeGqe^Pi^))' < (y^G(x,y)g,(y)G(y,z)g,(z)rfy)'rfz^02(z)rfz, 

by Cauchy Schwarz. This shows that 

\\QqeQqe\\\uj) < j^^ (y^G(x,y)g,(y)G(y,z)rfy)'g,2(z)dzdx. 

When ge(z, u) is bounded P— a.s., the above proof leading to (fT2l) apphes and we obtain 
(fT3!) under hypothesis [HI]. 

Using Lemma 12.11 we obtain that 

E{g.(y)ge(Og'(z)} < (^^)^^(O) + (^^)^^ (^T^) " 
Under hypothesis [H2], we thus obtain that 

mGQsQqef} < y^^G(x,y)G(x,C)(/^^(^^)G(y,z)G(C,z)dydCrfxdz 

Because r'^~^ip^{r) is integrable, then x i— > </)5(|x|) is integrable as well and the bound 
of the first term above under hypothesis [H2] is done as in f|T^ by replacing -R(x) by 
(/3 2 (|x|). As for the second term, it is bounded, using the Cauchy Schwarz inequality, by 

since x i— (y9(|x|) is integrable, D is bounded, and (jlj) holds. □ 

Applying the previous result to the process ge(x, cu) = q{^,uj), we obtain from the 
Chebyshev inequality that 

P(-;l|a».6*ll>rt<5^^f^<.^. (15) 

On the domain G Q of measure P(f2e) < e'^ where H^^e^g'ell > p, we modify the 
potential and set it to e.g. 0. We thus construct 

g.(x,-) = |Q (16) 

We have 

Lemma 2.4 The results obtained for qs{yi, uj) = q{-,uj) in Lemma lKB hold for qJyi. uj) 
constructed in ([16]). 
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Proof. For instance, 

E{||^g,6;/f} = EWM||6;g,6;/f} + Eto\f^.H||6;g,^/f} 
= E{xn\n4^)meGfr} < Eimom < e'^l 



The same proof holds for the second bound (fT3|) . □ 

We also need to assume that the oscillatory integrals studied in subsequent sections are 
not significantly modified when q{-,uj) is replaced by the new ^^(x, uj). We assume that 



[H3] 



limE 

£^0 



£2 V 



e 



1 /X 

£2 Vfc 



0. 



Note that such a condition is automatically satisfied when e~"''g(x, u) is bounded P-a.s 
for < a < i. 

With the modified potential, ([8]) admits a unique solution P-a.s. and we find that 



heWiuj) <\\gf\\ + \\gqegf\\ p- a.s., 

where || ■ || denotes L^{D) norm. Using the first result of Lemma 

n\\uef]<\\fr. 

Now we can address the behavior of the correctors. We define 

Mo = Qf\ 

the solution of the unperturbed problem. We find that 

(/ - gqeQqe){Ue - Uq) = -QqeGf + QQaQqeGf- 

Using the results of Lemma 12.21 we obtain that 



(17) 



we find that 



:i8) 



(19) 



(20) 



Lemma 2.5 Let be the solution to the heterogeneous problem ([T]) and Uq the solution 
to the corresponding homogenized problem. Then we have that 



-of}) 



2 < 



(21) 



Note that if we write = A^/ and uq = AqJ', with A,, and Aq the solution operators of 
the heterogeneous and homogeneous equations, respectively, then we have just shown 
that 

E{\\A-Aof}<e''. (22) 
Now QqeQqeiue — Uq) is bounded by e'^ in L^iVt] L'^{D)) by Cauchy-Schwarz: 



m'3qeQqe{u,-u,)\\}<(E{\\gq,gq,\\^}) (E{\\u, - u,f]) <e'. 
We need the following estimate: 
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Lemma 2.6 Under hypothesis [H2] of Lemma \2.^ we find that 

nwQqeGqsGfw'} < e'"'^ wff « s'^ur, (23) 

where rj is such that ^ ^ J y)'^^ uniformly bounded on D. 

Proof. By Cauchy Schwarz, 

\gq,gq,gf{^)\' < (^^G(x,y)g,(y)G(y,z)g,(z)G(z,t)dydz)'dt. 
So we want to estimate 

A = E{f G(x, y)G(x, C)qe{y)qe{C)G{y, z)G(C, ^)g.(z)g,(^)G(z, t)G(£, t)d[£Cyzxt]}. 
We now use ([9]) to obtain that A < Ai + A2 + A3, where 

A, = J^^ G(x, y)G(x, C)^^ (^^) G(y, z)G(C, ^)^^ (^^) G(z, t)G(|, t)d[|Cyzxt] , 
^2 = j^^ (y^^G'(x,y)G(y,z)y.^(^^^)G(z,t)rfyrfz)'rftrfx, 

A3 = G(x, y)G(^, t)G(x, C)G(z, t)^^ (^^) G(y, z)G(C, |)^^ (^^^) ^[ICyzxt]. 
Denote Fx,t(y;z) = ^(x, y)G'(y, z)G'(z, t). Then in the Fourier domain, we find that 

Here ^^{p) is the Fourier transform of x 1— > (/?^(|x|). Since (p^(ep) is bounded because 
r'^~^(y92(r) is integrable on R"*", we deduce that 

Ai < e'"' [ G\^,y)G\y,z)G\z,t)d^dydzdt < e"", 
Jd* 

using the integrabihty condition imposed on G(x, y). 

Using 2ab < a'^ + b'^ for (a, 6) = (G(x, y), G(x, 0) and (a, 6) = t), G(z, t)) 

successively, and integrating in t and x, we find that 

As < j^^ G{y, z)G(C, ^)^^ (^^) V^^ (^^) ^IyC^^] , 
thanks to dH). Now with (a, 6) = (G(y, z), ^)), we find that 

^3 < J^^ G' (y , z) v.^ (^^) (^^) ^[yCz|] ^ ^ 



2d 



since (/?2 is integrable and G is square integrable on the bounded domain D. 
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Let us now consider the contribution A2. We write the squared integral as a double 
integral over the variables (y, z, ^) and dealing with the integration in x and t using 
2ab < a"^ + b"^ as in the ^3 contribution, obtain that 



Using Holder's inequality, we obtain that 



A. < 



with p = 2 + ?7 and p' = since (p^ (^r)r'^ ^, whence ip'^ (^r)r'^ ^, is integrable. □ 
The above lemma applies to the stationary process ge(x, cu), and using the same proof 
as in Lemma [2. 4[ for the modified process ^^(x, cu) in (ITB]) . We have therefore obtained 
that 

E{\\u,-u + gq,gf\\}<e''. (24) 
For what follows, it is useful to recast the above result as: 

Proposition 2.7 Let g(x, cu) be constructed so that [H2]-[H3] holds and let ge(x, cj) be 
as defined in (1161) . Let be the solution to and Uq = Qf . We assume that uq is 
continuous on D. Then we have the following strong convergence result: 



limE 



Ue - Up 

d 
£2 



£2 



0. 



(25) 



Proof. Thanks to hypothesis [H3], we may replace ge(x, cj) by ge(x, cj) = q(^,uj) in 
2M up to a small error compared to e^. Indeed, 



—g(q\^-,ujJ-qs{-,uj)juo 

^ ^ 1 

< ||^||||Mo||L-'(D)E|xn=(c^^ 



52 Vfc 



< 1. 



■\gq(-,u)uo 

£2 



The rescaled corrector £~'iQq{-,uj)uQ does not converge strongly to its limit. Rather, 
it should be interpreted as a stochastic oscillatory integral whose limiting distribution 
is governed by the central limit theorem [TH [22] . We consider such limits first in the 
one-dimensional case and second for arbitrary space dimensions. 



2.2 Oscillatory integral in one space dimension 

In dimension li = 1, the leading term of the corrector e~'^{ue — uq) is thus given by: 



Uie{x,u) 



1 y 

-G{x, y)—=q{-, uj)uo{y)dy, 
D \Je e 



(26) 



where D is an interval (a, fe). The convergence is more precise in dimension d = \ than 
in higher space dimensions. For the Helmholtz equation, the Green function in c? = 1 is 
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Lipschitz continuous and we will assume this regularity for the rest of the section; see 
the next section for less regular Green's functions. Then Uis{x, u) is of class C{D) P-a.s. 
and we can seek convergence in that functional class. Since uq = Qf, it is continuous 
for / e L^iD). 

The variance of the random variable Uii;{x,uj) is given by 

E{m^^(x, cj)} = / G{x,y)G{x, z)-r(- ^uo{y)uo{z)dydz. (27) 

Jd'^ eve/ 

Because R{x) is assumed to be integrable, the above integral converges, as e ^ 0, to 
the following limit: 

E{ul{x,uj)}= [ G\x,y)R{0)ul{y)dy, (28) 
Jd 

where 

R{r)dr = 2 E{q{0)q{r)}dr. (29) 

-OO Jo 

Because (l26l) is an average of random variables decorrelating sufficiently fast, we 
expect a central limit-type result to show that Uis{x, u) converges to a Gaussian random 
variable. Combined with the variance (1291) . we expect the limit to be the following 
stochastic integral: 

ui{x,uj) = -a G{x,y)uo{y)dWy{uj), (30) 



D 



where dWyi^uj) is standard white noise on {C{D), B{C{D)),F) [S]. More precisely, we 
show the following result: 

Theorem 2.8 Let us assume that G{x,y) is Lipschitz continuous. Then, under the 
conditions of Proposition 2.7. the process Uis{x,uj) converges weakly and in distribution 
in the space of continuous paths C{D) to the limit Ui{x,uj) in f l30p . As a consequence, 
the corrector to homogenization satisfies that 

£ dist. ^ _^ / (^(^2;, ?/)Mo(?/)(iVFy, as e — > 0, (31) 



D 



in the space of integrable paths L^{D). 

Proof. We recall the classical result on the weak convergence of random variables 
with values in the space of continuous paths [S] : 

Proposition 2.9 Suppose [Zn] 1 < < 00) are random variables with values in the 
space of continuous functions C{D). Then Zn converges weakly (in distribution) to Z^ 
provided that: 

(a) any finite-dimensional joint distribution (Z„(xi), . . . , Zn,{xk)) converges to the joint 
distribution (Zoo(xi), . . . , Zoo{xk)) as n ^ 00. 

(b) {Zn) is a tight sequence of random variables. A sufficient condition for tightness 
of {Zn) is the following Kolmogorov criterion: there exist positive constants v, 13, 
and 5 such that 

{i) supE{|Z„(t)|''} < 00, forsometeD, 

n>i (32) 
(n) E{\Zn{s) - Zn{t)n < \t - s\'-^^ 

uniformly in n > 1 and t,s E D. 
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Tightness. Tightness of Uis{x,u!) is obtained with u = j3 = 2 and 6=1. Indeed, we 
easily obtain that 

E{|ni,(x,a;)p} <1, 
in fact uniformly in x G -D. Now by assumption on G{x,y) we obtain that 

E{\uUx, u) - uU^, = e( y) - G{^, y)]l=q{^)uo{y)dy^ ' 

= [ [G{x,y)-G{^,y)][G{xX)-G{^,0]-R{^)My)uo{C)dydC 



< |x - / -\R{^)\u,{y)uo{C)dydC <\x- 
Jd2 e e 



since the correlation function R[r) is integrable and Uq is bounded. This proves tightness 
of the sequence Uie{x, u), or equivalently weak convergence of the measures generated 
by on {C{D),B{C{D))). 

Finite dimensional distributions. Now any finite-dimensional distribution {ui^{xj, ^))i<j<r 

has the characteristic function 

$,(k) = ^{e^k,uu{x,,u^)^^ k= {ki,...,kn). 

The above characteristic function can be recast as 

n 

As a consequence, convergence of the finite dimensional distributions will be proved if 
we can show convergence of: 

/ m{y)-^q{-)dy Im-= [ m{y)adWy, e ^ 0, (33) 

J D ^ JD 



'■me 



for arbitrary continuous moments m{y). Such integrals have been extensively analyzed 
in the literature, see e.g. [9l [32], where the above integral, for D = {a,b) may be 
seen as the solution Xe{b) of the following ordinary differential equation with random 
coefficients: ^ 

Xe = —=q{-)m{t), Xe{a) = 0. 
A/e e 

Since we will use the same methodology in higher space dimensions, we give a short 
proof of (1331) using the central limit theorem for correlated discrete random variables as 
stated e.g. in [TO] . 



Approximation by piecewise constant integrand. Note that if we replace m{y) 
by mh{y), then 

E{(/„,,-/„,,)2} < ||m-m,||L, (34) 

where || ■ ||oo is the uniform norm on D. It is therefore sufficient to consider (133!) for a 
sequence of functions rrih converging to m in the uniform sense. Since m is (uniformly) 
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continuous, we can approximate it by piecewise constant functions rrih that are constant 
on M intervals of size h = Let rrihj be the value of m/j on the j**^ interval and define 



the random variables 



Mej = TUhj 



1 y 

{3-l)h ^ 



Independence of random variables. We want to show that the variables M^j be- 
come independent in the limit e ^ 0. This is done by showing that 



M 



^(k) = E{e* ^^=1 ''^^^^ } - JJ Eje^'^^*^^^ } 



as e — > 0, 



for all k = {kj]i<j<M e M^. Let k e M*^ fixed, < ?7 < | and define 

J{j-l)h+ri ^ 

Now we write 

Using the strong mixing condition ((61), we find that 



Now we find that E{Q^^} = and E{[Q^^]^} < 77. The latter result comes from in- 
tegrating e"^R{^)dsdt over a cube of size 0{rf). Since |e*^' — 1| ^ we deduce 
that 



for an arbitrary random variable Z (equal to or to X]fi2 ^j^ej here). Thus 



By induction, we thus find that for all < 77 < I, 



, kiM, 



This expression tends to say for t] = 62. This shows that the random variables M^j 
become independent as e — * 0. We show below that each M^j converges to a centered 
Gaussian variable as e — *• 0. The sum over j thus yields in the limit a centered Gaussian 
variable with variance the sum of the M individual variances. 
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Central Limit Theorem for discrete random variables. By stationarity of the 
process we are thus led to showing that 

/ -^qi-)dy [ adWy = aWh = (J^^{0, h), e ^ 0, 

Jo ^ Jo 

where Af{0, h) is the centered Gaussian variable with variance h. We break up h into 
N = h/e (which we assume is an integer) intervals and call 

2 y ri 
Qj = / -(l{-)dy = / q{y)dy, j G Z. 

J(j-l)e ^ ^ Jj-1 

The Qj are stationary mixing random variables and we are interested in the limit 

j=l V 

Following remark 3 in [10], we introduce Am and as the a— algebras generated by 
{lj)j<m and {qj)j>m, respectively. Let then 

f ^ f E{(r/-E{r/})(e-E{e})} ^ ,2, , , . , .2, ,nx\ .o.^ 
p{n) = sup <^ ^ 1 r] e L [Ao), L [A }\. (36) 

Then provided that ^^nyiPi^) ^ obtain the following central limit theorem 

rr N 

7= J2 ^ VhaAfiO, 1) = aXiO, h), (37) 

where Af{0, 1) is the standard normal variable, where = is used to mean equality in 
distribution, and where = Xlnez ^{lo^ln}- K remains to verify that the two definitions 
of a above and in (12^ agree and that J2n>iP(^) ^ Note that 

^E{gogn} = / / E{g(t/)g(z)}dyrfz= / / E{q{y)q{y+z)}dydz = / ^(0)d|/ = ^(0), 

Jo Jo J-00 Jo 

thanks to (!29|) . Now we observe that p(n) < (f{n — 1) so that summability of p(n) is 
implied by the integrability of (/?(r) on M+. This concludes the proof of the convergence 
in distribution of Ui^ in the space of continuous paths C{D). 

It now remains to recall the convergence result (1251) to obtain (13 ip in the space of 
integrable paths. □ 

2.3 Oscillatory integral in arbitrary space dimensions 

In dimension 1 < d < 3, the leading term in the corrector e~^(Me — Uq) is given by: 

uie(x,a;)= I -G{:>^,y)^qs{y,uJ)uo{y)dy. (38) 
15 



The variance of uie(x, tu) is given by 

E{ul^{x,Lj)} = j ^G{x,y)G{x,z)^R(^^^—^yuoiy)uo{z)dydz. 

As in the one-dimensional case, it converges as e ^ to the hmit 

E{ulix,co)} = a' [ G'ix,y)uliy)dy, a' = [ E{q{0)q{y)}dy. (39) 
Jd Jr'' 

Because of the singularities of the Green's function (^(x, y) in dimension d > 2, we 
prove here less accurate results than those obtained in dimension c? = 1 in the preceding 
section. 

We want to obtain convergence of the above corrector in distribution on [Q, JF, P) 
and weakly in D. More precisely, let Mfc(x) for 1 < k < K he sufficiently smooth 
functions such that 

mfc(y) = - / M,,(x)G'(x, y)uo{y)dx = -gMk{y)uo{y), 1 < k < K, (40) 
Jd 

are continuous functions (we thus assume that uo(x) is continuous as well). Let us 
introduce the random variables 

4e(^) = / mk{y)\q(-,uj)dy. (41) 
Jd 62 ^£ ^ 

Because of hypothesis [H3], the accumulation points of the integrals hei^) are not 
modified if q{-, u) is replaced by ^^(y, u). The main result of this section is the following: 



Theorem 2.10 Under the above conditions and the hypotheses of Proposition \277[ the 
random variables Ike{'^) converge in distribution to the mean zero Gaussian random 
variables Ik{uj) as e — >■ 0, where the correlation matrix is given by 

S,fc = E{J,-4} = (T^ / m,(y)mfc(y)dy, (42) 
Jd 

where a is given by 

a'= [ E{q{0)q{y)}dy. (43) 

jRd. 

Moreover, we have the stochastic representation 



hiiu) = / mkiy)adWy, (44) 
'd 



where dWy is standard multi-parameter Wiener process I3^ . 
As a result, for M(x) sufficiently smooth, we obtain that 



(^^^,m) ^-a [ gM{y)gf{y)dWy. (45) 
^ £2 / Jd 
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Proof. The convergence in (H5|) is a direct consequence of (jH]) since 



D2 



M(x)G(x,y)Mo(y)c^WyC^x= / gM{y)gf{y)dWy, 

Jd 



and of the strong convergence fl25l) in Proposition 12 .71 The equahty (jUj) is directly 
deduced from fH2]) since /^(ci;) is a (multivariate) Gaussian variable. In order to prove 
( 1421) . we use a methodology similar to that in the proof of Theorem 12. 8[ 
The characteristic function of the random variables Ikei'^) is given by 

$,(k) = E{e^S"-i'=^^-(-)}, k = (A^i, . . . , kK), 
and may be recast as 



-d 



K 



$^(k) = E{e^^o™(^)"^^(^'")'^y}, m(y) = ^kjmj{y). 

i=i 

So (142|) follows from showing that 

h{uj)= [ miy)^q(^,uj)dy^ [ m{y)adWy, (46) 

for an arbitrary continuous function m(y). As in the one-dimensional case and for the 
same reasons, we replace m(y) by mh{y), which is constant on small hyper-cubes Cj of 
size h (and volume h"^) and that there are M ^ h~'^ of them. Because dD is assumed to 
be sufficiently smooth, it can be covered by Ms ~ h~'^~^^ cubes and we set mh(x) = 
on those cubes. The contribution to lei^^) is seen to converge to as /i — > in the 
mean-square sense as in flM|) . 

We define the random variables 

Mej{uj) = mhj I 4r?(-, ^)dy, 1 < J < M, 
Jc, 8 2 e 

where ruhj is the value of rrih on Cj and are interested in the limiting distribution as 
£ ^ of the random variable 

M 

/»= J]M,,(^). (47) 
j=i 

We show below that these random variables are again independent in the limit e ^ 
and each variable converges to a centered Gaussian variable. As a consequence, I^{uj) 
converges in distribution to a centered Gaussian variable whose variance is the sum of 
the variances of the variables M^j{u)) in the limit e — > 0. 

That the random variables M^j are independent in the limit e — > is shown using a 
similar method to that of the one- dimensional case. We want to obtain that 



M 



S{k) = E{e' ^^=1 ''^^^^ } - Y[ Eie"'^^^^ } ^ as e ^ 0, for all k = {kj}j e 
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Let < r/ < I and V] = {x G Cj; ci(x, dCj) > r]}. We define 
We write again: 

Using the strong mixing condition iQ, we find that 

We find as in the one-dimensional case that E.{Q^j} = and E{[Q^^]^} < rfh^'^"^^ < 77 
with a bound independent of e. This comes from integrating e~'^R{^^—^)d'x.dy on a 
domain of size 0{[irih'^~^Y). The rest of the proof follows as in the one-dimensional case. 

It remains to address the convergence of M^j as e — > 0. By invariance of q'(x), it is 
sufficient to consider integrals on the cube [0, h], with h = (/i, . . . , /i). It now remains 
to show that 

/ ^q(l,uj)dy^a [ dWy = aX{0, h''). (48) 

^[0,h] £2 V£ / J[o,h] 

For a multi-index j 6 Z*^, we define 



/j+[o,i] 

Then (148|) will follow by homogeneity if we can show that 

^ E (49) 

^'^^ je[o,n] 

The latter result is proved in e.g. [lOl [19]. The results in these references are stated 
in terms of a-mixing coefficients. Since a coefficients are bounded by p coefficients [HI 
p. 4], we state the results in terms of less optimal p- mixing coefficients. 

Let A and B be subsets of Z'^ and let A and B be the a algebras generated by gj on 
A and B, respectively. Then we define 

(^ f E{(r/-EM)(e-E{0)} .^ryr^x MAm-> \ 

p[n) = sup < — J T] e L [A), ^ G L {B\, d{A, B) >n>. 

^ (E{r/2}E{e})^ ^ 

We then assume that E{g?} < 00 as in hypothesis [H2] and that p{n) = o{n~'^) and that 

00 

J^n'^-V^(n) < 00. (50) 

n=0 

Then we verify that the hypotheses in [10] (see also [19], p. 48]) are satisfied so that (H9l) 
holds with 

a2 = ^E{gogj}. 

We verify as in the one-dimensional case that the above a agrees with that in definition 
( 1431) . Now we verify that (!50|) is a consequence of the integrability of r'^~^(y9 2 (r). The 
decay p{n) = o{n~'^) is obtained when ^{r) decays faster than r~'^~^ for some 77 > 0; see 
Remark O □ 
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2.4 Larger fluctuations, random and periodic homogenization 

We now consider several generalizations of the results presented in earlier sections and 
compare homogenization in periodic and random media. The results stated in the pre- 
ceding sections, corresponding to the case a = below, generalize to larger fluctuations 
of the form: 

1 /X 

The corrector —Qq^Qf is now of order £"^(1"") for < a < |. The next-order corrector, 
given by QqeQqeGf in fl20l) . is bounded in L^{Q x D) by e'^*^2+^~^"^ according to Lemma 
12.61 The order of this term is smaller than the order of the leading corrector £"^(2-") 
again provided that < a < 2{2+ri) ' ^hich converges to | for (i=l,2as?7^oo and 
converges to | for d = 3 as r] ^ 1. 

In dimensions d = 1,2, we can infer from these results that e~'^^^"°'\ue — uq) 
converges in distribution to the limits obtained in the preceding sections as e — 
provided that < a < |. The proof presented in this paper extends to the val- 
ues < a < |. Indeed, the proof is based on imposing that the spectral radius of 



GQeGQe is sufficiently small using (1131) in Lemma 12.21 which for (ISTj) . translates into 
E{||^ge^ge||^^^2(£,))} ^ ^d{i-ia) ^ ^^^^ then verify that all results leading to Proposition 

12.71 generalize when < a < | to yield (!25|l with ei replaced by A proof of 

convergence for < a < ^ would presumably require us to analyze all the terms in the 
formal expansion 

00 

Ue = Y.^-QqefQf. (52) 

fc=0 

something we do not address here. In the limiting case a = |, the above theory breaks 
down and Ue no longer converges to the deterministic solution uq as is shown in the 
temporal one- dimensional case in |43j . 

The results on the corrector — uq obtained in the preceding sections, namely 
Theorems 12.81 and 12.101 are valid for 1 < d < 3. If we admit that the expansion fl52l) 
involves a first term Mq, a second term —GqeUo, and smaller order terms, then the 
results obtained in Theorem 12.101 show that — uq converges weakly in space and in 
distribution to a process of order 0(52). The critical case d = 4 yields a correction of 
order e^, whereas would be even smaller for d > 5. 

The theory presented in this paper does not allow us to justify (152|) when d > 4 
because the corresponding Green's function are no longer square integrable. Another 
argument shows that corrections of order correspond to a transition and that we 
should not expect quite the same results for c? < 3 and d > 4. Indeed, let us consider 
the problem in the periodic case: 

Ue = dD, 

on a smooth open, bounded, domain D C M'^, where g(y) is [0, 1] '^-periodic. Then 
following [7j, we introduce the fast scale y = - and introduce a function = U£(x, -). 
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Gradients Vx become ^Vy + Vx and fl53|) becomes formally 

( - ^Ay - ^Vx ■ Vy - Ax + g(y))M,(x, y) = /(x). 

Plugging the expansion = mq + eui + e'^U2 into the above equality and equating like 
powers of e yields three equations. The first equation shows that Uq = Uo(x). The 
second equation shows that Ui = Mi(x), which we can choose as Ui = 0. The third 
equation — Ay'U2 — AxMq + Q'(y)'Wo = /(x), admits a solution provided that 

-AxMo + (g)Mo = /(x), D 

with mq = on dD. Here, (q) is the average of q on [0, 1]'', which we assume is sufficiently 
large that the above equation admits a unique solution. We recast the above equation 
as uq = Qof ■ The corrector U2 thus solves 

-AyUs = {^{q) - g(y))no(x), 

and is uniquely defined along with the constraint (^2) = 0. We denote the solution 
operator of the above cell problem as so that U2 = —Q^{q — {q))Qf- Thus formally, 
we have obtained that 

n,(x) = ^/(x) (g))(^)^/(x) + l.o.t. (54) 

We thus observe that the corrector U2e(x) := n2(x, ^) is of order 0{e^) in the sense, 
say. In the sense of distributions, however, the corrector may be of order o(e™) for all 
integer m in the sense that M(x)u2£(x)(ix <^ for all m when M(x)uo(x) G C^{D). 
A similar behavior occurs for the random corrector 

vie{x,uj) = J -G{x,y)q(^^,u^uoiy)dy. (55) 

Theorem 12.101 shows that (fie,M(x)) is of order 0{e^) for M(x) and mo(x) sufficiently 
smooth and that e~2 (y^^, M(x)) converges in distribution to a Gaussian random vari- 
able. This result, however, does not hold in the L^(D)— sense for d > 4 when G(x, y) is 
the fundamental solution of the Helmholtz equation —A + go(x) on D. Indeed, we can 
prove that 

Proposition 2.11 Provided that mo(x) and -R(^) are sufficiently smooth, we obtain 
that: 



ER(x,^)} 



e''R{0) J G\x,y)uliy)dy l<d<3 
eVnel^^^^ulix) d^4 (56) 

•^4 



e%l{x){2i:Y f ^d^ d>5. 
Here ~ fc^ means = bs{l + o(l)). 
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Proof. We calculate: 

E{vl{^,uj)} = y^^G'(x,y)G(x,z)i?(^^)no(y)no(z)dydz. (57) 
Extending G(x, ■) by outside of D, by the Parseval equality this is equal to 

where J^x^^ is the Fourier transform from x to In dimension 1 < < 3, since R{e$,) — *■ 
R{0) pointwise, the Lebesgue dominated convergence theorem yields the result. In 
dimension d > A, however, the Green function is no longer integrable and the integral 
is larger than e"^. 

Let us consider the cases d> 4. We first replace G(x, y) by c^jx — y p"*^ where Cd is 
the measure of the unit sphere 5*^"^. The difference G(x, y) — qIx — yp"'' is a function 
bounded by C|x — y]^"*^, which yields a smaller contribution to E{i;f^}. We leave the 
details to the reader. We also replace Uo(y) by Mo(x), up to an error bounded by |x — y|" 
as soon as Uq{x.) is of class C°'"(-D). This contribution again provides a lower order term 
to K{vl^}. Similarly, we replace mo(z) by mo(x) and thus obtain that 

Let a > and -B(x, a) the ball of center x and radius a so that -B(x, a) C D. Because 
all singularities occur when y and z are in the vicinity of x, we use the proof of the case 
1 < (i < 3 to show that up to a term of order e'^, we can replace D by i?(x, a) so that 

EK(x,^)}-n^(x) / ^ ^,l_.^,i,. ^(^)^y^z- (58) 

JB^{0,a) Cd\y\ Cd\Z\ \ e / 

Now for (i > 5, using the dominated convergence theorem, we can replace -8(0, a) by 
because the Green function is square integrable at infinity, whence 

This, however, by the Parseval equality, is equal to 

EK,(x,..)}~n^(x)(2vr)'^ / ^e'^^O^^ = «o(x)(2vr)'^ / ±.e'mdt 

since the Fourier transform of the fundamental solution of the Laplacian is |^|~^. 

When d = 4, we come back to fl55]l . and replace one of the integrals on -8(0, a) by 
an integral on M'^ using again the dominated convergence theorem. This yields the term 

JB{0,a)xR^ C4|y| |Z| \ 8 / JB{0,a)xR'' C4|y| |^| 

= / ^f''£M ^)e'^-''d^dy = m{2iTer [ ^,dy + O(e^) 

JB{0,SL)xR'i C4|y| Isl Jb{0,^) C4|y| 
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Here, we have assumed that |-R(^) — -R(0)| was bounded by C|^|^ for some (3 > 0. □ 

In all dimensions, we thus obtain that £:~2fi£(x, u;) converges (weakly and in dis- 
tribution) to a limit ui(x, u;) = — JjjG{'x,y)uo{y)dWy. In dimensions 1 < d < 3, we 
have proved that ui was the limit of the corrector to homogenization e~^{ui; — Uq). The 
above calculation shows that the limit ui captures all the energy in the oscillations of 
the homogenization corrector vig in the sense that the limit of the norm E,{vl^{x, v)} 
is equal to the norm K{ul{x, v)}. 

In higher dimensions d > 4, as in the case of homogenization in periodic media, some 
energy is lost while passing to the (weak) limit. The corrector ui^ = e~2Vis converges 
weakly and in distribution to the limit ui. However, while the energy of the limiting 
corrector is £2 (E{||ui||^2(£)}2 , the energy of the true corrector (E{||fie|||2(£)}2 is of order 

0{e'^) for d > 5 and of order 0(£:^| lne|5) for d = 4. Most of the energy of the correctors 
is lost in passing from to its weak limit ui. 



3 Correctors for one- dimensional elliptic problems 

In this section, we consider the homogenization of the following one- dimensional elliptic 
problems: 

~'^^'^(^^^'>^^^ + (9o + qe{x,uj))ue = pe{x,u;)f{x), xeD = {0,l), ^^g^ 
uM = u,il) = 0. 

We consider homogeneous Dirichlet conditions to simplify the presentation. The coef- 
ficients as{x,u) and Pe(x, a;) are uniformly bounded from above and below: < < 
air(x,uj), Ps{x,u) < The (deterministic) absorption term go is assumed to be a 

non- negative constant. The generalization to a non- negative smooth function qQ{x) can 
be done. 

We assume that as{x, u) = a(|, u), qe{x, lj) = g(|, u), and Pe{x, u) = p(|, u), where 
a{x,uj), q{x,uj), and p{x,Lj) are strictly stationary processes on an abstract probability 
space {Q,J-',F). We will modify the mean-zero process ^^(x, w) as in the preceding 
section and assume here to simplify that q{x,u!) is bounded P— a.s. We also assume 
that the cross-correlations -R/g(x) are integrable for {/,(?} G {a,q,p}, where 

Rf^i^) = E{f{y,uj)g{y + ^,uj)}. (60) 

We also assume that the coefficients are jointly strongly mixing in the sense of ([6j), where 
for two Borel sets A and S in M*^, we denote by JF4 and jFg the a-algebras generated 
by the random fields a(x, a;), g(x, a;), and p(x, cu). We still assume that the p-mixing 
coefficient (p{r) is integrable and such that ip2 is also integrable. 

In the case where Qe = and = 0, the corrector to the homogenization limit uq 
has been considered in [12]. For general sufficiently mixing coefficients with positive 
variance cr^ = 2 K{a{0)a{t)}dt > 0, we obtain that Ue — uq is of order y/e and 
converges in distribution to a Gaussian process. This section aims at generalizing the 
result to ( l59l) using the results of the preceding section and a change of variables based 
on harmonic coordinates [M] . 
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Let us introduce the change of variables 



1 dze ^ a* 1 

a,(t) ' dx a,{x)' " E{a-i}' ^ ' 

and Ueiz) = ^^(x). Then we find, with x = x{zs) that 

-{a*)'^-^u, + a*qoU^ + ae[{l-a^^a*)qo + q^]u^ = a^pJ\ < z < z^{l) ^^^^ 

UeiO)=leiz,{l))=0. 

Let us introduce the following Green's function 

-a*^G'(x, y; L) + q^G^x, y;L) = 5{x- y) ^g^^ 



Then, defining 
we find that 



G(0,2/;L) = G'(L,y;L) = 0. 
qe{x,uj) = (1 - a~^{x,uj)a*)qo + ^^(x,^;), (64) 



Jo 

We recast the above equation as 

Ue{x,Uj) = Qeipef - QeUe), QeUix) = G{Zs{x) , Ze{y)] Zs{l))u{y)dy . (65) 

^0 

After one more iteration, we obtain the following integral equation: 

Ue = GePef " GeQeGePef + QeqsGsqeUe- (66) 

Since aoa*x < z^{x,uj) < a*aQ^x P— a.s., the Green's operator is bounded P— a.s. 
and the results of Lemma [2.21 generalize to the case where the operator replaces Q. 
As in ( IT5|) . we thus modify (i.e. we modify and q^) on a set of measure e so that 
ll^e^e^e^ell < T < 1 and assume that [H3] holds. 
Let us introduce the notation 

Ps = p + Sps, p = E{p}, ge = G + sge, g = E{g,}, uo = gpf. (67) 

We also define 

5z,{x) = z,{x)-x = / b(-)dt, b{t,uj) = ——--l. (68) 

Jo a{t,uj) 

We first obtain the 
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Lemma 3.1 We have that 

E{\6z,{x)6z,{y)\}<e, 0<x,y<l. (69) 
The operator may be decomposed as 

Ge = G + Gle + ne, (70) 

where 

/I r\ r\ 

{Sz,{x)— + S^e{y)g^ + Sz,{1)—)g{x, y- l)f{y)dy. (71) 

We also have the following estimates 

mGier}<e, E{||7^,||}<e. (72) 
Proof. We first use the fact tliat 

n\5z,{x)5zM\} < {^{{5z,{x)5zMf}Y ■ 
Denoting by he{x,uj) = b{^,uj), we liave to sliow tfiat 

Jo Jo Jo ^^^^l)^^^'^2)&e(2:3)fee(^4)c?[^l 2:2^3^4] } 

Now using tlie mixing property of the mean- zero field feg and the integrability of V9^(r), 
we obtain the result using (Q as in the proof of Lemma 12.61 

The integral defining is split into two contributions, according as ?/ < x or y > x. 
On these two intervals, G{x, y; L) is twice differentiable, and we thus have the expansion 

d d d 

G{z,{x), z.iy); z,{l)) = G{x, y; 1) + {s^eix)— + ^z^iy)-^ + 5z,{l)-^G{x, y; 1) + r„ 

where the Lagrange remainder = r^{x, Zi;{x),y, z^iy), z^{l)) is quadratic in the vari- 
ables {Sze{x), 6ze{y), 6ze{l)) and involves second-order derivatives of G{x, y; 1) at points 
(^,C,^) between {x,y;l) and {zs{x), Zs{y); Ze{l)). 

From (169!) and the fact that second-order derivatives of G are P— a.s. uniformly 
bounded on each interval y < x and y > x (we use here again the fact that aQa*x < 
Ze{x,u!) < a*aQ^x P— a.s.), we thus obtain that E{|r£(.)|} < e. This also shows the 
bound for E{||7?.e||} in (172]) . The bound for E{||^ij|p} is obtained similarly. □ 

Because we have assumed that and were bounded P— a.s., we can replace by 
Q -\- in (|66|) up to an error of order e in L^{Q; L'^{D)). The case of and bounded 
on average would require to address their correlation with defined in the proof of the 
preceding lemma. This is not considered here. 

We recast (1^ as 

Ue-Uo= {Gepe -Gp)f - GeQeGePef + GeQeGeqeiUe - Uq) + GeqsGeqeGf- (73) 

Because G^Zir^x), Zs{y); z^{l)) and p^ are uniformly bounded P— a.s., the proof of Lemma 
12.21 generalizes to give us that 

mGeqeGeqeW'} + E{\\GsqeGePsf\\'} + E{||(^,p, - Gp)ff} < S. (74) 

So far, since moreover ll^eO'e^eO'ell < r < 1, we have thus obtained the following result: 
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Lemma 3.2 Let Ue he the solution to the heterogeneous problem (159|) and Uq the solution 
to the corresponding homogenized problem. Then we have that 

{n\\ue-u4'})HV~s\\f\\. (75) 

The estimate ( 122|) with c? = 1 is thus verified in the context of the elhptic equation ( l59l) . 
As a consequence, we find that EjH-Ug — mqH^} ^ b so that by Cauchy Schwarz and (174|) . 



^{WQeqeGeQeiUe " Mo)||} < £■ 

It remains to exhibit the term of order ^/e in m^ — mq. Let us introduce the decomposition 

Ue-Uo= Qiep + Q5pe - QQeGp f + S^, (76) 
Se = {SGeSpe + T^sp) f " {QeqeQePe - QqeQp)f + QeqeQeqe{Ue - Uq) + QeqeQeqeQ f ■ 



Lemma 3.3 Let f G L^{D). We have 

E{b.||}<^||/||. (77) 

Proof. Because (^(^^(x), ^^(y); ^^(l)) is uniformly bounded P— a.s., the proof of 
Lemma \T6\ generahzes to show that IE{||^ege^ege^/P} < e^H/p. We aheady know 
that E{||7?.e||} < e. It remains to address the terms Ji = Qi^dpef, h = GqeQiePef, 
h = GieqeGpsf, and h = QqeQ^Pef- 

Because is uniformly bounded P— a.s., the first three terms are handled in a 
similar way. Let us consider E{/^}, which is bounded by a finite number (three here) 
of operators of the form 

E{ / 5ze{v^{x, y))H{x, y)6z,{v2{x, C))H{x, ()6ps{y)6ps{()f{y)f{C)dxdydC}, 

where Vk{x,y) is either x, y, or 1 for /c = 1,2, and H{x,y) is a uniformly bounded 
function. Using the definition of Sz^, we recast the above integral as 

/ / / E{b,{t,)be{t2)6pe{y)6p,{C)}dtidt2H{x,y)H{xX)f{y)f{OdxdydC. 

J D'i Jo Jo 

Using iQ , we see that the above integral is bounded by terms of the form 

rJo Jo V'K^^^T^)^K^^^7^)^^^^^^'^^'^'^^^^'^'^^''^^^^''^^^^'^^^^^ 

where (mi, ti2, ^^3, ^4) = {ui,U2,U3,U4)(ti,t2,y,C) is an arbitrary (fixed) permutation of 
iti,t2,y,C)- Because (/?(r) is integrable, the Cauchy Schwarz inequality shows that the 
above term is < The term E{J|} is given by 

e{ J^^ g{x, y)G{x, CMy)qe{OG{y, z)G{C, OSps{z)Spe{Of{z)fm[xyza]]. 
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Since G{x, y) is uniformly bounded on D, we again use (Q as above to obtain a bound 
of the form 11/ Ip. □ 

It remains to analyze the convergence of the contribution [QuP + Q^Pe ~ G(jeSp]f- 
As in fl2BD. we define 



We recast the above term as 



1 



GleP + G^Pe - GqeQp fix). 



b[-)Hh{x,t) + 6p[-)Hp{x,t) -q[-)Hg{x,t) dt, (79) 



with 



Hp{x,t) 
Hg{x,t) 



Xx{t)—G{x, y; 1) + Xyit)—G{x, y; 1) + -7^G{x, y; 1) pf{y)dy 



dy 



dL 



G{x,t)f{t) 

G{x,t) [ G{t,z)f{z)dz. 



(80) 



where Xx{t) = 1 H < t < x and vanishes otherwise. We have the following result. 

Theorem 3.4 Let f G L°°{0, 1). The process Uii,{x,u!) converges weakly and in distri- 
bution in the space of continuous paths C{D) to the limit ui{x,u!) given by 



Ui{x,uj) = / a{x,t)dWt, 
Jo 

where Wt is standard Brownian motion and 

POO 

a2(x,t) = 2 / E{F{x,t,0)F{x,t,T)}dT, 
F{x, t, t) = Hb{x, t)b{T) + Hp{x, t)6p{T) - H,{x, t)q{T) 
As a consequence, the corrector to homogenization thus satisfies that: 

Ue ~ Uq dist. . X „ 

^ — [x) >ui{x,uj), as £ — > U, 

in the space of integrable paths L}{D). 
We may recast Uie{x,uj) as 

3 



^11 



^2) 



^3) 



Uie{x,u) 



/ p,i-)Hkix,t)dt, 
k=i e 



where the pk are mean-zero processes and the kernels Hk{x,t) are given in (IHOjl . The 
corrector in (jSTl) may then be rewritten as 



3 

^lix) = / (rk{x,t)dWi 



^5) 
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with three correlated standard Brownian motions such that 

dW^dWl" = pjkdt, (86) 

where we have defined 



Pjk 



J 


poo 

1 Efe(0)pfe(r) 




+ Pk{0)Pj{T)}dT 




poo 

E{p,iO)p,iT)}dT j 


PCO 1 

1^ E{pk{o)Pk{T)}dT}y 



17) 



That ( IHTi) and ( l85l) are equivalent comes from the straightforward calculation that both 
processes are mean zero Gaussian processes that have the same correlation function. 
The new equation (!85|) shows more clearly the linearity of the ak, whence Ui{x), with 
respect to the source term f{x). 
Proof. We recast Uis{x,uj) as 



Uis{x,Uj) 



/ qk{-)Hk{x,t)dt, 



with a different decomposition as in ( l84l) . where the qk are mean-zero processes and the 
kernels Hk{x,t) are given implicitly in (IHOl) . We verify that we can choose the terms 
Hk{x,t) in the above decomposition so that all of them are uniformly (in t) Lipschitz 
in X, except for one term, say Hi{x,t), which is of the form 

d 

Hi{x,t) = Xx{t)Li{x,t), Li{x,t) = j —G{x,y;l)pf(y)dy, 

where Li{x, t) is uniformly (in t) Lipschitz in x. This results from the fact that G{x, y; 1) 
is Lipschitz continuous and that its partial derivatives are bounded and piecewise Lips- 
chitz continuous; we leave the tedious details to the reader. 

Because of the presence of the term Hi{x, t) in the above expression, it is not sufficient 
to consider second-order moments of uu as in the proof of Thm. 12.81 Rather, we consider 
fourth-order moments as follows: 

E{Mx,u;)-uU^,u)n= 4 E / ^{qkA-)Qk/-)qk/-)qk/-)}y< 

£ J D"^ £ £ £ £ 

fci,fc2,fc3,fc4 

4 

n {Hk^{x,tjn) - Hk^{^,tm))dtidt2dt3dU. 

m=l 

Using the mixing condition of the processes qk and Lemma 12.11 (where each g in ([9]) 
may be replaced by qk without any change in the result), we obtain that E{\uis{x, u) — 
uislC,, u;)|^} is bounded by a sum of terms of the form 



4 

r it t it t' 

/ if^- -)ipH- -)Yl(^kmix,tm) - HkM,tm))dtidt2dt3dt4, 



m=l 
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whence is bounded by terms of the form 

When all the kernels if^^ are Lipschitz continuous, then the above term is of order 
\x — ^1^. The largest contribution is obtained when k\ = k2 = \ because Hi{x,t) is not 
uniformly Lipschitz continuous. We now concentrate on that contribution. We recast 

Again, the largest contribution to the fourth moment of comes from the term {Xx{t) — 
X^{t))Li{x,t) since Li{x,t) is Lipschitz continuous. Assuming that x > without loss 
of generality, we calculate that 

iXx{t) - xdi))Liix,t){xds) - xds))Lii^,s)-ipk- — -)dtds 

1)2 EE 

-X fX 



I I Li{x,t)Li{^,s)-ipk^-^)dtds<{x-0, 
Je Je EE 



since is integrable. Note that this term is not of order — xp. Nonetheless, we have 
shown that 

E{\m,{x,u)-m,i^,u)\^}<\^-x\^, 

so that we can apply the Kolmogorov criterion in Prop. 12.91 with u = 2, j3 = 4, and 
6 = 1. This concludes the proof of tightness of uis{x, ui) as a process with values in the 
space of continuous functions C{D). 

It remains to verify step (a) of Prop. 12.91 The finite-dimensional distributions are 
treated as in the proof of Thm. I2.8l and are replaced by the analysis of random integrals 
of the form: 



-f 

Jo 



dt. 



The functions m are continuous and can be approximated by rrih constant on intervals 
of size h so that we end up with M independent (in the limit e ^ 0) variables of the 
form: 

nibhbj + mphdpj + nighqj. 



N . ^ 



It remains to apply the central limit theorem as in the proof of Thm. 12.81 The above 
random variable converges in distribution to 



Af{0, ha ), a = 2 E{ {nibhb + rriphSp + mqhq) (0) {rribhh + mph5p + mquq) {t)}dt. 
Jo 

This concludes our analysis of the convergence in distribution of to its limit in the 
space of continuous paths C{D). The convergence of — uq follows from the bound 

dZZD. □ 



28 



4 Correctors for spectral problems 



4.1 Abstract convergence result 

For uj ^ Q, let Aj^{u) be a sequence of bounded (uniformly in u P— a.s. and in 77 > 0), 
compact, self-adjoint operators, converging to a deterministic, compact, self-adjoint, 
operator A as 77 —> in the sense that the following error estimate holds: 

E|| A(cj) - A\\P < r]P, for some 1 < p < 00, (88) 

where ||A^(ct;) — A\\ is the L'^{D) norm and D is an open subset of W^. 

The operators A and P— a.s. A^(co') admit the spectral decompositions (A„,m„) and 
(A^,M^), where the real- valued eigenvalues are ordered in decreasing values of their 
absolute values and counted m„ times, where rrin is their multiplicity. 

For A„, let /in be (one of) the closest eigenvalue of A that is different from A„. Let 
us then define the distance: 

dn = (89) 

Following [31], we analyze the spectrum of A^j in the vicinity of A„. Let F be the circle 
of center A^ and radius dn in the complex plane and let R{(,A) = {A — be the 
resolvent of A defined for all complex numbers ( ^ cr(^); the spectrum of A. The 
projection operator onto the spectral components of B inside the curve F is defined by 



^n[fi] = -— J^R{C,B)dC. (90) 

Note that for all ( E T, we have that R{(, A)Pn[A] = (A^ — C)~^- We can then prove 
the following result: 

Proposition 4.1 Let A^ and A be the operators described above and let A„ be fixed. 
Then, for 77 sufficiently small, there are exactly rrin eigenvalues A^ of A^ inside the 
circle F. Moreover, we have the following estimates: 

n\K - + E{|||< - unV) ^ J A 1, (91) 

for a suitable labeling of the eigenvectors of A^ associated to the eigenvalues A^. 

Proof. It follows from [HH Theorem IV. 3. 18] that for those realizations u such that 
II A^(co') — A\\ < dn, then there are exactly m„ eigenvalues of in the (i^— vicinity of F. 
Since this also holds for every Am such that dm > dn, we can index the eigenvalues of 
Ari as the eigenvalues of A. Moreover, 

\Xl{u;)-Xn\<\\A,{u)-A\\. 

For those realizations u such that ||A^(co') —A\\ > dn, we choose m„ eigenvalues of A^(co') 
arbitrarily among the eigenvalues that have not been chosen in the (im— vicinity of Am 
for |Am| > |A„|. 

For all realizations, we thus obtain that 

\X1{UJ) - Xn\ < - 



dn 
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It remains to take the pth. power and average the above expression to obtain the first 
inequahty of the proposition. 

In order for the eigenvectors and m„ to be close, we need to restrict the size of t] 
further. To make sure the eigenvectors are sufficiently close, we need to ensure that 

- P4A] = ^ Jm. A) - ^(C, A)]dC = ^ ^ ^(C, A,){A, - A)R{C, A)dC 

is sufficiently small. On the circle F and for ||y4 — y4^|| < dn, we verify that 

sup||i?(C,A)|| = j-, sup||i?(C,A,)|| < ] 

by construction of dn and by using R~^{(, A^) = -R~^(C, A) + (A;, — A) and the triangle 
inequality 

\\R-\C,A,)\\ > \\R-\CA)\\-\\A,-A\\>dn-\\Ar,-A\\. 
Upon integrating the expression for P„[y4^] — P„[A] on F, we find that 

p := \\P4A,] - P4A] II < ■ L"^,?;!'' < l-ll A - ^il < 1, 
for 2\\An - A\\ < d^. 

For self-adjoint operators A and A^, the above bound on the distance p between the 
eigenspaces is sufficient to characterize the distance between the corresponding eigen- 
vectors. We follow [311 1.4.6 & II. 4. 2] and construct the unitary operator 

Wn ={l- {Pn[Ar^ - Pn[A]fy~' (P„K]Pn[A] + (/ - Pn[A,])(J - P„[A])). (92) 

Let Un,ki ^ ^ k < rrin be all the eigenvectors associated to an eigenvalue A„, n > 1. 
Then the eigenspace associated to admits for an orthonormal basis the eigenvectors 
defined by [21] 

K,k = Uyn,k, l<k<mn. (93) 
The relation fl92l) may be recast as 

U:i = {I- Rl){l + P„[A,](P„[A,] - P„[A]) + (P„[A,] - P4A])P„[A,]), 

where ||P^|| ^ p^. This shows that 

ll^n-^II^P and \\ul^k-Un,k\\^P^^\\Av- A\\, I < k < m^, 

whenever (i~^||A^ ^ ^11 < for sufficiently small. When (i~-^||A^ ~ ^11 > Z^? "we find 
that \\un,k — ^^fcll ^ 2yu||A^(u;) — A||/(i„, where the vectors are constructed as an 
arbitrary orthonormal basis of the eigenspace associated to A^. Upon taking pth power 
and ensemble averaging, we obtain fl9T]) . □ 
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4.2 Correctors for eigenvalues and eigenvectors 

Let {Xn,Un) be a solution of Aun = XnUn and let and be the solution of A^m^ = 
A^M^ defined in Proposition 14.11 We assume that flSSj) holds with p = 2. 
We calculate that 

^^^^ = ^^Y^^n) + - ^n, ((A, - A:1) -{A- A„))n„). 

The last term, which we denote by is bounded by 0{ri) in L^{Q) using the results 

of Proposition HIT] with p = 2 and the Cauchy Schwarz inequality. Thus, r'!^{uj) converges 
to in probability. 

Let us assume that the eigenvectors are defined on a domain D cM.'^ and that for a 
smooth function M(x), we have: 

(M(x), ^''~'^ n„(x)) / M(x)a„,(x,y)rflVyrfx as r/ ^ 0. (94) 

Using this result, and provided that the eigenvectors Un(x) are sufficiently smooth, we 
obtain that 



Xn ~ Xn dist. ^ 



/ Un{x)aniyi,y)dWydx:= / A„(y)dWy as r/ 0. (95) 
Jd2 Jd 



The eigenvalue correctors are therefore Gaussian variables, which may conveniently be 
written as a stochastic integral that is quadratic in the eigenvectors since (T„(x, y) is a 
linear functional of Un- The correlations between different correctors may also obviously 
be obtained as 

Xn ~ A„ — Am "1 V^O 



E -11 ^ ^ / A„(x)A„(x)dx. (96) 

Let us now turn to the corrector for the eigenvectors. Note that 

\\Un - = 2(1 - {Un,ul)), 

SO that {un,Un) IS equal to 1 plus an error term of order 0{rf) on average. The con- 
struction of the eigenvectors in (jHSD show that n„ — u1 is of order 0{rf) in the whole 
eigenspace associated to the eigenvalue A„. It thus remains to analyze the convergence 
properties of (m„ — m^, Um) for all m ^ n. A straightforward calculation similar to the 
one obtained for the eigenvalue corrector shows that 

(^^^^, {A - A = - ( (A-a:0-(a-a„) ^^^ ^ _ 1 _ ^^^^^ _ ^^^^ 

The last term converges to in probability (and is in fact of order 0{r]) in L^{Q) as 
above). We thus find that 

f<^^um) ^ ^ [ Um{^)an{x,y)dWydx. (97) 

The Fourier coefficients of the eigenvector correctors converge to Gaussian random 
variables. As in the case of eigenvalues, it is straightforward to estimate the cross- 
correlations of the Fourier coefficients corresponding to (possibly) different eigenvectors. 
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4.3 Applications to some specific problems 



The first application pertains to the following problem: 



Lemma 12.51 and its corollary (l22l) show that (IHSjl holds with p = 2 and rj = . The 
operators A,, and A are also compact and self-adjoint for a large class of operators 
P(x, D) which includes the Helmholtz operator P(x, D) = —A + qo{x.). 

Let (A^,M^) be the solutions of X^P^u^ = and (A„,m„) the solutions of XPu = u. 
Then we find that 



d 
52 



A„a / Un{^)G{^,y)u^{y)dWydi^=-\la I ul{y)dWy, (99) 



or equivalently, that for the eigenvalues of and P, we have: 



£2 



a / ul{y)dWy. (100) 



The Fourier coefficients of the eigenvectors satisfy similar expressions. 

The second example is the one-dimensional elliptic equation (l59l) . Still setting 77 = 
£2^ we find that 

" Un / o-„(x,t)rfWt, 

where o"„(x, t) is defined in (182|) with the source term / in (IHOj) being replaced by m„(x). 
The operators A^ and A satisfy (IHH!) with p = 2 thanks to Lemma 13.21 and its corollary 
(I22!) . The expressions for the eigenvalue and eigenvector correctors are thus directly 
given by (l95l) and ( 1971) . respectively. 



4.4 Correctors for time dependent problems 

As an application of the preceding theory, let us now consider an evolution problem of 
the form 

Ut + ePu = 0, t>0, m(0) = mo, (101) 

where e is a constant, typically e = 1 or e = and P is a symmetric pseudodifferential 
operator with domain T>{P) C L'^{D) for some subset D d and with a compact 
inverse A = P~^, which we assume without loss of generality, has positive eigenvalues. 
We then consider the randomly perturbed problem 

+ eP^Ur, = 0, t > 0, M^(0) = Mo, (102) 

where Pri{oj) verifies the same hypotheses as P with compact inverse = P~^. 

We assume that A^ and A are sufficiently close so that flSSl) holds. Following the 
notation of the preceding section, we denote by A„ and the eigenvalues of A and 
and by u„ and u^^ the corresponding eigenvectors. 

We then verify that 

u{t) = e-'^^uo = e-'^"\un, %)«« := Yl "«(^)^«' ""W = e"'^"*(^^n, %)• 
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and 

n 

We can now compare the Fourier coefficients as follows: 



{un, uo) + e-'^-\^ no) + r„ (103) 

where |r^| ^0 strongly in Lp{Q) as ?7 — »■ 0. This may be recast as 

= e-^^"*6t^^^(n„,no) + e-^^"*(^^^^, no) + s,, (104) 

T) 7] 7] 

where — > strongly in Lp{Q) as r/ ^ 0. 

The above difference thus converges to a mean zero Gaussian random variable whose 
variance may easily be estimated from the results obtained in the preceding section. 

Since we do not control the convergence of the eigenvectors for arbitrary values of 
n (because we do not control in this study the speed of convergence in distribution 
of the random correctors), we cannot obtain the law of the full corrector n,,(t) — u(t). 
We can, however, obtain a corrector for the low frequency parts M^(t) and u^(t) of 
Uri{t) and u{t), respectively, where only the N ffist terms are kept in the sum in the 
index n. We may easily estimate the corrector for {u}^(t) — up^(t),Um) using the above 
expansion for the Fourier coefficients and the results obtained in the preceding section. 
We again obtain that the corrector is a mean zero Gaussian variable whose variance 
may be calculated explicitly. 

Other time-dependent equations may be treated in a similar way. For instance, the 
wave equation 

utt + Pu = 0, n(0) = no, Ut{0) = go, (105) 

where P is a symmetric operator with compact and positive definite inverse, may be 
recast as 

Wt-Aw = 0, w{0)=wo, ), ^=(o (106) 



utj ' 

We verify that the eigenvalues A„ of A are purely imaginary and equal to ±zVXp, where 
Xp are the positive eigenvalues of P. The orthogonal projector onto the nth eigenspace 
of A is found to be 

■np,_A2 
Anp_A2 



so that 



P-A2 W^O 

\Up.2J [go 



(107) 



A similar expression may be used for the perturbed problem Urjit) , where P is replaced by 
Prj. The results presented earlier in this section easily generalize to provide an estimate 
for the low frequency component of u(t) — Ur,(t). We leave the details to the reader. 

When the Green's function associated to the operator dt + eP is sufficiently regular, 
for instance when eP = —A, more refined results may be obtained by considering 
expansions similar to the expansion ([8]) considered for steady-state problems. We do 
not consider such developments here. 
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5 Conclusions 



We have considered the corrector to the homogenization of source and spectral problems 
for the Helmholtz equation with highly oscillatory random potential. The method works 
because the operator Qq^ appearing in ([8]) may be seen as lower-order, in the sense that 
it converges rapidly to with e. This requires that the homogenized solution ^/ be a 
good approximation to the source problem ([8]). The method was then generalized to 
the one-dimensional elliptic problem flS^ . which after a change of variables to harmonic 
coordinates, may also be recast as an integral equation fIBB]) with a term ^^g^ that may 
also be seen as lower-order. 

Such expansions are not currently available for more challenging problems of the 
form —V ■ as{'x,uj)'Vus = /, augmented with appropriate boundary conditions. The 
use of the Green's function to the homogenized elliptic equation does not allow for a 
rapidly converging expansion of the form or fl66p . The analysis of correctors for such 
equations, for which current state of the art estimations are given in [IS], remains an 
open problem; see [17] for a related discretized elliptic equation. 

The correctors were analyzed here in the setting where the random coefficients have 
integrable correlation function -R(x) in ([5]) (and additional mixing properties). The ex- 
pansions in ([H]) and fl^Bl) may be generalized to random coefficients with correlation func- 
tions -R(x) which decay as Ixl"""^ for some < a < 1 as |x| ^ oo. In such frameworks, 
following the expansions obtained in [5], we expect random correctors with Gaussian 
statistics and amplitudes of order rather than , at least for dimensions 1 < < 3 
for the Helmholtz problem. These long-range effects will be analyzed elsewhere. 
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